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INTRODUCTION 


The effect of finite deformation upon a rotating, orthotropic cylinder 
was investigated by Sandman (ref. 1). He was able to predict the influence of 
finite deformations and relate his results to the degree of orthotropy. In 
this study an attempt has been made to study the same problem using a general 
incremental theory. 

The incxremental equations of motion are developed using the variational 
principle discussed by Washizu (ref. 2). A more than adequate development of 
the goveraing equations has been given by Atluri (ref. 3). Although his inten- 
tion is to implement a finite element scheme to solve boundary value problems, 
the equations are given in general tensor notation. Hofmeister, Greenbaum, and 
Evensen (ref. 4) have presented an excellent discussion of the use of an incre- 
mental analysis; again, their goal is the application of a finite element anal- 
ysis. The governing equations are also developed in the treatise by Biot (bef . 
5), using both a geometrical viewpoint and a variational method. The governing 
equations are rederived here, in somewhat less detail, using the principle of 
virtual work for a body with initial stress (ref. 2). 

The governing equations are reduced to those for the title problem and a 
numerical solution is obtained using finite difference approximations. Since 
the problem is defined in terms of one independent space coordinate, the finite 
difference grid can be modified as the incremental deformation occurs without 
serious numerical, difficulties. The nonlinear problem is solved incrementally 
by totaling a series of linear solutions. This method was used to solve the 
same problem discussed in ref. 1 and gave identical results. 


GOVERNING EQUATIONS 


The derivation of the governing equations is based upon an incremental 
variational principle (ref. 2). The body is assumed to be in equilibrium at 
some arbitrary reference state along the load path. Let 

X = a + u (1) 


be the transformation of a particle at point a to point x in the same space, 
then u is_^the displacement of the particle. At the beginning of some increment 
of load, a is the initial coordinate and x is the current coordinate, and the 
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two are identical. Let initial stresses o^, initial surface tractions t°, and 
initial body forces act on the body before the addition of the load incre- 
ment. These stresses emd loads are with respect to the initial coordinate axis 
and are referred to a unit area before the loading increment is applied; hence, 
they are referred to an undeformed area and volume. 

Assuming the initial stress system is in equilibrium, it follows that. 


div a® + f ® = 0 

(2) 


0®n = 

(3) 


where n is a unit normal vector . If the body is then loaded with some increment 
of surface traction or body force, the total stresses at the end of that incre- 
ment of load are the sum of the initial stresses and incremental stresses. 

In order to formulate the principle of virtual work, first define a non- 
linear strain tensor, such as. 


where 


D = 

E + N 

(4) 

E = 

(Vu + Vu^) 

(5) 

N = 

(ViJvu) 

(6) 


where u is the displacement field corresponding to ^ and D_ is referred to as 
Green’s strain tensor (ref. 6). The notation is basically the direct notation 
used by Gurtin (ref. 7), although some symbols are different. 

Introduce a virtual displacement 5u and incremental stresses, body forces, 
and surface tractions, £, ?, and respectively. The principle of virtual work 
for a body with initial stress may be written. 


^{(o^+£)*6^-(f°+f)*Su}dv-| (t°+t)*6udS = 0 


(7) 


where Sq corresponds to the surface on which stresses are specified. Substitut- 
ing equations (5) and (6) into (7) and noting that ^ and o_ are symmetric yields 


f {o® • 6Vu+a° • Vu'^6 Vu+a • 6 Vu+a ♦ Vr^SVuJdv- /(f^ • 6u+f • 5u)dv- 
y— — V 

I (t°*6u+t*6u)dS = 0 

Making use of 18(1) (ref. 7), equation (8) can be rewritten as 
JSu'Cdiv a + div(a° Vu"^) + div(^Vii^) + ?]dv 


( 8 ) 


-/ 5u*[on + (a®Vu^)^n + (oVu^)^n - ?]dS = 
J 6u*[a°n - t°]dS - /6u*[div o° + $°]dv 


TsT^ 


(9) 
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According to equations (2) and (3) the righthand side of equation (9) should be 
zero; therefore, the equations of equilibrium become 

div u + div[(£° + q)Vu'^] + $ = 0 (10) 

and the boundary condition is 

Op + [(of + o_)Vu^3'^n = t (11) 

The assumption that the incremental strains are small Implies that u is small 
incrementally and 

£ = E, i.e. N = 0 in equation (4). (12) 

The initial stress may not be small; hence, we retain of terms in equations , (10 ) 
and (11). It follows from equation (12) that for a linear incremental stress- 
strain relation the incremental stress will be small. Therefore products of 
can be neglected and the governing equations become 

div a_ + div(£fVj^) + ^ = 0 (13) 

cn + = t (14) 

Equations (2) and (3) serve as an error check and can be used at any increment 
to determine the equilibrium status of the initial stress system. 


The total stress u_ at the end of any load increment becomes the initial 
stress of for the next load increment. Then, £ must be referred to the initial 
coordinates a and the deformed area ;Ln order to become of . The transformation 
has been given by Fung (ref. 6) and can be rewritten as 

o° - (p/p )V xoV (15) 

— o a-- a-^ 

where p/pq is the ratio of final mass to initial mass and indicates that the 
operator is with respect to the initial coordinates It follows from equation 
(1) that 


V X = V ( a+u ) = 6 +Vu (16) 

a—~ a a~ 

where is a unit tensor . For an incremental theory equation (16) may be written 
i + V^u = £ + Vu = J (17) 

It follows that 

p/p^ = det|Va| - 1 - tr(7u) (18) 

where tr( ) represents the trace of a tensor. Combining equations (15) through 
(18) gives the transformation 

£° = [1 - tr(7u)]J0j'^ (19) 
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where 


= (6 + Vu)"^ = 6 + Vu'’’ 


GOVERNING EQUATIONS FOR A ROTATING CYLINDER 


( 20 ) 


The general equations can be reduced to plane cylindrical coordinates in 
order to implement the analysis of a rotating cylinder. The problem is one of 
axisymmetric plane strain; hence, the displacement vector ^ reduces to Up, the 
radial component, which will be referred to as u. 

The numerical method will be applied to the equation of equilibrium (13), 
which in plane cylindrical coordinates may be written 

a' + (a - o^)/r + u' + a^(u” + uVr) ~ ag u/r^ t 
r r 0 r r 0 

pca^(r + u) = 0 (21) 

where f = p(r+u)m^ the inertia force, cr^ and oq are the radial and tangential 
stresses, respectively, and the prime denotes differentiation with respect to r. 

Equation (12) is represented by the linear strains 

E = u’ and En = u/r (22) 

r y 

Following Sandman (ref. 1) we assume a linear anisotropic stress-strain relation 

0 ^ = u’ + C ^2 (23) 

0g = 0^2 + C^2 (24) 

Substituting equations (23) and (24) into equation (21) yields the incremental 
governing equation 

u" + u’/r - au/r^ + o° + u’(o°* + o°/r)/C^j^ 

- og u/Cj^^r2 + po)2(r+u)/Cj^^ =0 (25) 

where 

a = ^ 22^^11 ^ " ^ 12^^11 

The boundary condition, equation (14), becomes 

u’(l + 0 °/C, J + 3 u/r = 0 (27) 

r 11 

The linearized incremental stress transformation, equation (19), becomes 

0 ° = 0 (1 + u' - u/r) (28) 

r r 
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= 0g(l - u' + u/r) 


(29) 


NUMERICAL ANALYSIS 


The governing equation (25) was solved using a finite difference technique. 
The primary constraint to be dealt with is the magnitude of each increment of 
strain. It must be small enough to insure that equation (12) is not violated. 
After each increment of displacement is calculated, the finite difference grid 
must be updated ; hence , the finite difference equations must be reformulated 
after each incremental solution. The difference operations may be derived as 
follows 


(du/dr). -(u.^j_-u.)/Ar2-(u.-u._^)/Arj_-(u.^^^ 


2u. 


2u. 

1 


(d^U/dr^) . --;p — y- r — ^ ^ ■ 

1 Ar,(ArT+Ai>_) Ar, Ar 


2’^ltl 


1-2 


(30) 

(31) 


The first incremental solution is merely the linear solution for the first 
increment of body force. Before the second incremental solution is determined, 
the initial stresses are assumed to be equal to the stresses obtained for the 
first increment. These stresses are transformed according to equations (28) and 
(29). The incremental displacement associated with each finite difference node 
is added to the coordinate of that node-; hence, a new initial stress problem is 
formulated. The nonlinear analysis for the equation developed by Sandman (ref. 
1) was obtained by transposing all nonlinear terms to the right. The displace- 
ments for the previous analysis were used to evaluate the naniinear terms, and 
a solution for u is obtained. The calculated displacements are then used to 
calculate new nonlinear terms, and the solution is repeated. This process con- 
tinues until the two sets of displacements agree to within some tolerance. 

This method was used to verify the results obtained by Sandman (ref. 1) and 
appears to be accurate and efficient. 


Equations (2) and (3) can be used at any Increment to determine if the 
initial stress system is still In equilibrium. If the initial stress system is 
not in equilibrium, the solution can be corrected by including equation (2) in 
the governing equation (25). 


NUMERICAL RESULTS 


Solutions were obtained for three different materials. These material 
parameters were assumed to approximate the behavior of steel, aluminum, and a 
composite epoxy- fiber orthotropic material. The maximum radial and tangential 
stresses are shown in figiire 1 as a function of The cylinder was assumed 

to have an outside radius of 0.127 m (5 inches) and inside radius of 0.254 m 
(10 inches). The maximum radial stress occurs approximately halfway between 
the inside and outside, while Cg is maximum at the inside radius. 
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The percent deviation of the nonlinear solutions •ahoye the linear is illus- 
trated in figures 2 and 3. The increase in stress using the equations of refer- 
ence 1 appear to be almost linear in every case. The radial stress increase, 
using the incremental theory, is similar for both steel and aluminum and • 
reflects a nonlinear behavior. The increase for the composite appears to become 
constant. The nonlinear tangential stress deviation increases and then tends to 
decrease for both isotropic materials; however, this behavior is not demon- 
strated for the composite. 

In all cases the increase in stress level does not appear to be significant 
for stresses in the elastic range. The analysis presented herein should be 
extended to include nonlinear material behavior. 
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